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Abstract 

The relativistic kinetic theory of the phonon gas in superfluids is developed. The technique of the 
derivation of macroscopic balance equations from microscopic equations of motion for individual 
particles is applied to an ensemble of quasi-particles. The necessary expressions are constructed in 
terms of a Hamilton function of a (quasi-)particle. A phonon contribution into superfluid dynamic 
parameters is obtained from energy-momentum balance equations for the phonon gas together 
with the conservation law for superfluids as a whole. Relations between dynamic flows being in 
agreement with results of relativistic hydrodynamic consideration are found. Based on the kinetic 
approach a problem of relativistic variation of the speed of sound under phonon influence at low 
temperature is solved. 
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I. INTRODUCTION 



Over the last years the problems of relativistic superfluid dynamics were discussed re- 
peatedly Such an interest is well justified. Superfluids turned out to be the object that has 
revealed some problems of relativistic hydrodynamics and thermodynamics in a new fashion. 



T 
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ed to earnest discussion in a relativistic generalization of Landau two-fluid model 

Usually, two approaches are used in relativistic hydrodynamics, Eckart and Landau- 
Lifschits ones. In the former the macroscopic velocity of a gas (or liquid) is defined as the 
unit vector parallel to the particle number flow, and in the latter it does as eigenvector of 
the energy-momentum tensor. In equilibrium all flows such as a particle number flow and an 
entropy flow are parallel to this unit vector. When dissipation effects are taken into account 
all the macroscopic values get additions that are orthogonal to the vector of the macroscopic 
velocity. 

Superfluids fundamentally differs from the perfect fluid. The particle number flow n % and 
the entropy flow S l as well as conjugated dynamic /ij and thermal Bj momenta have differ- 
ent directions even in equilibrium. It allows to develop an approach that does not prefer 
any direction. In context of this approach it becomes possible not only to describe ade- 
quately relativistic superfluid dynamics but also to discuss some problems of the traditional 
approaches 13,^1. 

nn 

According to p, |y] the energy-momentum tensor of superfluids at non-zero temperature 
can be represented in the form 

7? = n'fij + S% -P5}, (1) 
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does not demonstrate the symmetry 
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where P is the liquid pressure. Though the notation 
of the tensor in an explicit form, it really takes place y, 

Among the four flows, ri 1 , S\[i l , and G l , only two, n l and S l say, are independent. The 
other two are expressed as 



// = Brf + AS\ & = An 1 + CS\ (2) 

with three scalar function A, B, and C named the anomaly, bulk, and caloric coefficients 
respectively. Consideration of superfluidity supposes the coefficient A should be negative 
and the coefficients B and C should be positive p. 



Since there is no a preferred flow in the theory developed in fl, fl] the vector of the 
"hydrodynamic" velocity may be associated with any of them. The choice is made by 
convenience reason. For example, the velocity may be parallel to the particle number flow 
when a gas of identical particles is considered, whereas it may be direct along entropy flow 
if thermodynamic properties are studied. In this work there also exists a preferred vector 
V 1 . It has a sense of a superfluid velocity. In context of this paper the superfluid velocity 
plays a role of a velocity of the carrying medium where quasi-particles, phonons and rotons, 
are located. 

Relativistic superfluid hydrodynamics can be used also for solving more practical prob- 
lems in different areas apart from studying fundamental matter. This formalism is found 
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to be efficient to describe processes in neutron stars 
cosmological models 
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111 ] and to construct the 



Hydrodynamic approach is essentially macroscopic and therefore it does not take into 
account microscopic processes that have influence on the dynamic and thermal parameters 
(such as entropy, pressure etc.). It is also clear that hydrodynamics can not describe all the 
phenomena in superfluids. For example, usage of hydrodynamic equations is restricted by 
the condition lot <C 1, where u is the sound frequency, r is the free path time. This condition 
is not fulfilled at temperatures close to zero because free path time of the quasi-particles 
increases. In this temperature range system behavior is described by a kinetic equation. 

Meanwhile kinetic phenomena in relativistic superfluids remain beyond the scope of study. 
The superfluid kinetic theory concerns the problems related to the excitation gas, which sig- 
nificantly determines values of fluid parameters at non-zero temperature. When dissipation 
is discounted, the quasi-particle gas determines thermal effects in the first place. In hydro- 
dynamics influence of the quasi-particles is masked by macroscopic flows such as an entropy 
flow. In essence the unit vector in the direction of the entropy flow is considered as a 
macroscopic velocity of the quasi-particle gas. 

The kinetic theory of the quasi-particle gas for non-relativistic superfluids was developed 
in To generalize this theory for the relativistic case is the primary purpose of this 

study. Covariant kinetic theory of the quasi-particles naturally develops and supplements 
relativistic kinetic theory for the real mass and massless particles. 

One of the main goals of kinetic theory is to connect macroscopic variables with intrinsic 
microscopic parameters. With respect to superfluids this means that a quasi-particle contri- 
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bution into the superfluid dynamic parameters should be obtained. At the low temperatures 
we can take into account only phonons (i.e. excitations with the linear dependence between 
the energy and the momentum) as influence of other quasi-particles may be neglected. 

When we deal with quasi-particles in special or general relativity, a main difficulty is to 
define the macroscopic variables adequately. The usual definitions |l6[ are inapplicable here, 
since the quasi-particles being the massless move slower than light. In Sec. IHII this problem 
is solved using the Hamilton function for the quasi-particles. Structure of the function 
is based on the acoustic metrics introduced in Sec. |Hj This permits to modify the usual 
constructions of the energy-momentum tensor, the particle number flow and the entropy 
flow. In Sec. I1VI balance equations for the quasi-particle gas are derived from the kinetic 
equations. Simultaneous consideration of the phonon balance equation and conservation 
laws for all liquid allows to find a phonon contribution to the chemical potential and the 
pressure of the liquid. In Sec. IS the problem of a sound speed variation under influence of 
the phonon gas at low temperature is investigated for the Minkowski and Robertson- Walker 
metrics. 



II. ACOUSTIC METRICS 

In the classical Landau's theory [l^, Q a superfluid state is governed by quasi-particles. 
Each quasi-particle has a certain momentum and energy that are functionally related. Such 
a functional dependence is said to be an energy spectrum of superfluids. In the spectrum of 
liquid helium one recognizes two kinds of quasi-particles, named phonons and rotons that 
contribute maximally to thermodynamic properties. Phonons correspond to the linear part 
of the spectrum whereas rotons belong to the nonlinear part of the spectrum. Processes 
being due to rotons are inessential at a temperature close to zero and their contribution is 
neglected. Thus at the low temperature a quasi-particle gas may be considered as a purely 
phonon gas. 

To develop our theory it is convenient to introduce so-called acoustic or phonon metrics 
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= 9 ij + [jj - I J V*V*, (3) 
where V % is the superfluid velocity normalized by ViV 1 = 1, and the sound speed c„ at zero 
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temperature is defined by the relation 



Cp _ dP _ n|| dp 



(4) 



c 2 dp no dn 

in which n\\ = n'Vi, n* is the liquid particle number flow, po is the chemic al p otential, Pq is 
the pressure, and p is the energy density of the liquid at zero temperature j2l|. The acoustic 
metrics is constructed so that phonon momenta satisfy the relation 

t j PiPi = 0, (5) 

i.e. they are tangent vectors to the characteristic hypersurface of sound propagation in the 
medium. Along with the tensor 7f J the metric tensor with subscripts is involved 

,2 

_p 

c 2 



9ij = 9ij + ( % ~ 1 ) ViVj, (6) 



which satisfy the condition ^g^ = 5\. 

Notice that the acoustic metrics (j3J) can be formally derived as the special case of a 
broader class of a tensor corresponding to an arbitrary energy spectrum. Actually, we can 
define the tensor 



G" = </'■> -\ ( ^ - 1 ) r r\ 



where e = piV 1 is the energy, and 7r 2 = —A'^piPj = —(g^ — V % V^)piPj is the square of 3- 
momentum of a quasi-particle in the local rest frame of the medium. In this case G^PiPj = 0. 
If we will establish a functional dependence between e and 7r we obtain the covariant form of 
the energy spectrum. However in general, G l i can not be used as a metrics, since it depends 
not on only coordinates but also on momenta. Only for phonons with the condition 

e = ^vr, (7) 
c 

G l i can be considered as a metrics. 

The superfluid velocity appearing in the acoustic metrics sets the direction for the "su- 
perfluid momentum" pj = pVj, where p is the chemical potential of superfluids 0. The 
momentum pj is gradient of a scalar function (which is interpreted as a phase of a wave 
function in a quantum description of superfluidity): Pj = V,-0. This permits to derive the 
equation 

V m] = 0, (8) 



that demonstrates irrotational nature of superfmid moving. Contraction of (JBJ) with V 1 gives 
the relation 

liDVj = A}V,/i, (9) 

(D = V l Vi is the convective derivative) that can be interpreted as generalization of the non- 
relativistic equation for the superfluid velocity. The irrotational condition can be obtained 
in a more traditional form by convolution of (jHJ) with the projector A\ and using ©: 

A i[fc V^-] = 0. (10) 



III. MACROSCOPIC VARIABLES 



To describe the phonon gas it is necessary to introduce the usual set of macroscopic 
parameters, such as particle and entropy flow four- vectors and an energy-momentum tensor, 
which are constructed using a phonon distribution function f(x,p). If to follow the usual 
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way [lfil 1191 ] the macroscopic variables should be defined as 

iV' = 2 C /p7(x.rt 9(P °gr' ) ^, (Ha) 

^'^'''' y ft <ub) 

but these expressions are not applicable in this case since they do not carry required physical 
content. To illustrate this, let us consider the vector N % in the flat space-time. 
The component N° must have meaning of a number particle density, i.e. 

N° = cj f(x,p)d 3 p, (12) 

where the integration is to be taken over usual (three-dimensional) momentum space and 
Po is determined by (j5J). N a (a = 1, 2, 3) are the components of a particle flow three- vector. 
Therefore in the medium rest frame the phonon particle flow should be defined as 

N a = c p J ^-f(x,p)d 3 p. (13) 

5-function provides transformation of (jllaj) to the usual momentum space 

N* = c[ Q^f(x,p)d 3 p. (14) 
J T Pk 
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Comparing this expression with (fT2j). one can conclude that N % defined according (jllaj) can 
not be used particle number flow. 

To give a physical sense to the macroscopic currents, their definitions should be modified. 
For this purpose we introduce Hamilton function for a quasi-particle 

n = \thiPi, (is) 

that is equal to zero at the hypersurface defined by the expression (jSJ). Using this function 
we construct macroscopic variables 

# = c t ™ t{x , r) *>^#L, (16a) 

J dpi {2-nh) 6 y/^g 



i dU p , ,e(p 7 -V J 'Wft) d 4 p 



S' = -ck B /^(/'"/ 



( 1 + /)M 1 + />]«^|. 

Further we will use a compact notation 

dV = e{e)5{H) d 4 p/{27ch) 3 ^j. 

The macroscopic flows defined in (|T6|) satisfy all the wished requirements. To verify this 
fact let us consider the modified vector iV\ Integration of (jlfiaj) with 5-function leads to the 
expression (JH|), but now the acoustic metrics Tf- 7 (not g^) is in the numerator, providing 
fulfilment both (IT2J) and dU). 

The determinations (fTTj) usually include 9(p°) that gives the positivity of the particle 
energy. When we deal with quasi-particles, a role of the energy is assigned to the convolution 
PiV 1 = e. That is why it fill the place of the argument to the 0-function in the expressions 

(HE). 

The expressions (116)1 can be considered as generalization of the traditional definitions of 
the macroscopic flows. To revert to the usual notations, Hamilton function should be taken 
in the form 

W = \ W 3 PiPj -m 2 c 2 } 



for massive particles, and 



W = ^g lJ PiPj 



for massless ones. 
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The entropy flow definition may be different according to particle statistics 
in superfluids are bosons, that leads to the construction in the square brackets of (|16c|) . 

The energy-momentum tensor defined by (j!6b|) is not symmetric in the gene ral case. 
This implies that the angular momentum of the phonon gas is not conserved |23j|. This is 
a quite expected circumstance because the gas is not an isolated system. The liquid is not 
only a medium for phonons, its parameters are influenced by the quasi-particles. In view 
of this interrelation, the energy-momentum of the phonon gas is also not conserved as we 
shall see later. Conservation laws hold for the system as a whole. It will allow us to make 
conclusion about a phonon contribution in the thermodynamic parameters of superfluids. 
The angular momentum for all the system is implied to be conserved, therefore the total 
energy-momentum tensor can be represented in the symmetric form (| . 

The macroscopic flows ()16|) can be naturally split into contributions parallel and orthogo- 
nal to V 1 . It is convenient to decompose the phonon energy- momentum tensor as following: 

T] = wV l Vj + Q% + Q 3 V l + IT} + QjV\ (17) 

where 

w = T}V i Vi= C - I e 2 f(x,p)dV, 
Qi = T 3 k V j A ki = c J eiCif(x,p) dV, 

II} = If A*AJ = C y n^f(x,p)dV, 
e = PiV\ TTi = A*p k . 

A similar decomposition one can perform for the entropy flow: 

5 i = s N y i + t7 i , s\\ = S i V i , a % = A)S j . (18) 

IV. PHONON CONTRIBUTION INTO SUPERFLUID DYNAMICS 
A. Kinetic equation for phonons 



The phonon distribution function satisfy Liouville equation 

df df OH df OH 



24, 
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ds dx l dpi dpi dx % 



J(x,f), (19) 
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where canonical Hamilton equations are used 

dx i OH d Pi OH 



(20) 



ds dpi ds dx l 

and J(x, f) is a collision integral for the phonons. 

The derivative of the Hamilton function (|15|) with respect to coordinates can be converted: 



an i dg kl - 

5* = 2 PkPl ^ = ~ PkPs9 1 



~ VmVk dp- m (21) 
Then the equation (|19|) takes the form: 

V = ™% f= J lx , p) , (22, 

where 

v. = ^+rLp^ (23) 

is the Cartan derivative |25J] and are the Christoffel symbols defined with respect to the 
metrics 'g 13 . 

One can directly verify that the following equalities are valid: 



as well as 



V iPk = 0, (24) 

V*^ = 0, (25) 
dp k 



1 dyf—g 1 dy/^g i 1 <9c p 



<7 ax* ^/ — g ox 1 c p ox 1 



Taking into account these properties it is easy to show that 

f 9n ~ , . — r m r , , ^ i dc p r an r , . 

J W P^y k f^p)*P = V k J W Pn--P,f^p)*P--— k J — Pil .. Pi J(x,p)dV, 

(27) 

where Vjt is covariant derivative defined with respect to the metrics "g^. 
B. Balance equations of phonon gas 

The relation (|27j) is very convenient to derive macroscopic equations. Let us multiply the 
Liouville equation (}2~2*|) on Pi 1 ..pi l and integrate it over momenta. If number of the momenta 



9 



I = then we obtain the balance equation 

ViJV* = / J{x,p)dP. (28) 



Notice that the covariant derivative is associated with the metrics g y and a rate of change 
of the quasi-particles is governed by the collision term of the kinetic equation. Therefore 
the equation (J28|) becomes a conservation law only in equilibrium when J(x,p) = 0. 

There is a different situation for phonon energy-momentum balance, that can be obtained 
setting 1 = 1. In this case the integral of the right hand side of the Liouville equation turns 
to zero by a momentum conservation law and hence we find that 

v /z] + (r™ - r™)^ = o, (29) 

where and are the Christoffel symbols with respect to the metrics g^ and g^ cor- 
respondingly. The second term in (}2"9*j) is found using the explicit expression for and 
decomposition ((T7j). Then the balance equation (}2T?j) can be written as 

(r 2 \ • w 
- - 1 QWjVi - -V jCp = 0. (30) 
c p J c p 

This means that there is no a phonon energy-momentum conservation even in equilibrium. 
Contracting the equation (}3T)j) with and applying the decomposition (fT7j) . we find 

w 

Dw + w9 + ViQ 1 - Q l DVi - U^ViVj Dc p = 0, (31) 

Cp 

where the following notations are introduced: 

9 = ViV* = 5% = A%, % = A fc(i V%. (32) 

Next we follow Carter 0, who noticed that there is no evidence of parallelism of entropy 
and particle flows in the hydrodynamic equations even in equilibrium. Formally we may 
construct hydrodynamics with unparallel equilibrium flows. Really it is realized in superfluid 
hydrodynamics. In usual liquid the particle and entropy flows become unparallel in non- 
equilibrium only, but traditionally a basic thermodynamic relations (that are related to a 
rest frame of the velocity selected) remain the same. Such an approach seems to be not 
well consistent since the flows have different directions, and hence, one more independent 
parameter appears in addition to a usual set of variables. It is an absolute value of difference 
of the flows. Hence the Gibbs relation have to be supplemented by a term of a kind XSY, 
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where Y is a new variable and X is a thermodynamically conjugated current. Using a 
Legendre type transformation the new parameters may be added in the Gibbs relation as 



— Y5X. If Y = yiy 1 then the variation of the scalar can be transformed to the variation 
of the vector: X5Y = (Xy i /Y)5y' t = Xi8y\ 

In our consideration the superfluid velocity and entropy flow have different directions and 
the Gibbs relation will be used in the following way 

5w = T5s\\ + cr l 5ri + Hi5n\\, (33) 

where T is the temperature, 73 is the current conjugated to a 1 . The third term appears 
because the phonon energy depends on the liquid particle density. The parameter /xi is not 
the chemical potential of the phonon gas (it equals to zero since the number of quasi-particles 
is not conserved) and it should be determined as well as r». 

Taking into account ()18|) the entropy production density can be realized in the form 

V i S i = Ds\\ + S||0 + Via 1 . (34) 

The liquid particle number flow can be decomposed at components parallel and orthogonal 
to the superfluid velocity 

n* = n\\\r + H\ n|| = nVi, H* = A)n j . (35) 

The corresponding balance equation may be written as 

Virf = Dn\\ +n\\6 + ViW = 0. (36) 

Dw is substituted in the equation (pfTjl using (|33|) . Expressing the convective derivatives 
Dsn and Dnu from ()34)1 and we find 



TV i S i -Ts\\9-TV i a i + aWT i -fi 1 n ll 9~fi 1 V i H i +w9 + V i Q i -Q l DV i -W j V l V J -—Dc p = 0. 

Cp 

(37) 

The speed of sound c p is assumed to depend only on n\\, and Q l = Qa\ where the 
proportionality factor G is to be defined. On simple rearrangement, the expression (|3T|) can 
be written as 

TV*? = |lF - A ij (w - T S| , - m] + -^g^J + j % + 

{T - 0} V,,^ + /ii - V^ 1 + {£{QVi - n) - Vi0} a\ (38) 
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In the expression (|38|) we have turned from the convective derivatives to the Lie derivatives 
with respect to the vector V\ This is due to the fact that the convective derivation is 
concerned with changes of vector projections onto direction V\ whereas the Lie derivative 
does changes of the vectors themselves. The derivatives of two types are related by the law 

CVi = DV h Cn = Dr % + t%j. (39) 

Each term on the right-hand side of (|38j) includes two factors. One of them consists of 
space-time derivatives, whereas another algebraically depends on the macroscopic variables. 

In equilibrium the entropy production density is equal to zero. Since spacelike components 
of the flows S l and n l are independent, each term should be put to be zero separately. This 
leads to 

IF = A ij (w - Ts [{ ) - oV, (40) 
T = 6, (41) 

ft = (42) 
c p on\\ 

v,e - c(QVi - n ) = o. (43) 

The first three conditions give the connection between variables involved in the energy- 
momentum tensor (j!7)l . The last one is non-dissipative heat conduction equation. 

In non-equilibrium the entropy production density is positive due to additional contribu- 
tions in the currents. It is usually supposed that for small deviations from equilibrium these 
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contributions depend only linearly on the gradients of the equilibrium quantities 
According this approach the right hand side of the equation (}4*3*|) becomes proportional to 
the non-equilibrium component in a 1 , corresponding to an irreversible heat flow. Since this 
dependence emerges in the next order of magnitude with respect to the equilibrium values, 
it does not affect on the relations (j4*mi - (|4*2^) (we disregard viscous effects). 

In the context of the discussion regarding infinite propagation velocity for heat signal 
in relativistic thermodynamics, Carter suggested [6| another way to obtain the positivity 
of the entropy production. This approach is based on the assumption that heat inertia is 
nonzero and the energy-momentum tensor is constructed as a purely algebraic function of 
independent currents such as a particle number flow and a heat current. This statement looks 
like Fok's "physical principle" which implies that energy-momentum tensor components must 
contain the state functions only and may not include gradients as well as coordinates in the 
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explicit form |26j. Taking into account Carter's ansatz, the relations (|40j) - (|%2j) remain valid 
as before and the equation (}4"3*|) transforms into 



where is a positive semidefinite resistivity matrix. The equation (}4"4"|) leads to a hyperbolic 
equation for heat conduction automatically without using an additional phenomenological 
parameter. 

For our purposes the result regarding relations (|40j ) -([42 p is more necessary than differ- 
ent forms of the heat conduction equation and it is the same at least for the considered 
approaches, however the later seems to be more consecutive in our case. 

C. Phonon contribution into superfluid dynamics 

In previous section phonon gas was considered. If it is in solids, then it is believed 
that phonon gas does not affect the moving of the medium. In superfluids, phonons not 
only determine thermal parameters but also change dynamic parameters of the medium. 
Thus considering superfluids as closed system containing two parts, liquid (medium) and 
phonons, we should assume that tensor (JT7jl is contained in the total energy-momentum 
tensor (JTJ). Because of is the energy- momentum tensor of all the system it obeys the 
covariant conservation law 




(44) 




(45) 



77 = n* H - P f 8) 



(46) 



where Pf is the pressure of the liquid properly, without a phonon influence. 
Using (|4T)|) the energy-momentum tensor ()17j) can be represented in the form 



T] = S l Q* - P ph 5). 



(47) 



In this case 



Tc 2 



(48) 




(49) 



P ph = Ts|| - w. 



(50) 
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Since the thermodynamic properties of superfluids are completely determined by the 
phonons, and the entropy flow vector appears in the second term of (0) only, it is reasonable 
to identify O* with Qj. 

The residual part has the form (J46)) with P p h + Pf = P, but parameters appearing in 
(I46|) change their values under the phonon influence. To obtain the phonon contribution 
into the dynamic parameters one should eliminate the phonon balance equation (|5U|) from 
the conservation law ([45)1 . that yields 

(c 2 \ w 
--l)Q i V j V i + -V j c p = 0. (51) 
c p J c p 

One transforms this expression using the condition ©: 

"i*H£)-*( p '-?fe) + 

fxW - - Q^j % = 0- ( 52 ) 

The most obvious way to provide the fulfillment of this relation is to equate the last term 
to zero separately This means that vectors H l are proportional to a 1 : 

litt { = - lj Ta\ (53) 

In the residual part the expressions in parenthesis should be identified as the chemical 
potential /io and the liquid pressure Pq at zero temperature that are connected by nyJ/xo = 
SPo- 

Such an approach is not only possible but it allows the equation (J52j) to be carried out for 
an arbitrary state equation. Thus, at the temperature close to zero the chemical potential 
and the pressure of the superfluid are determined by the equations 

These expressions coincide exactly with the nonrelativistic ones obtained for the phonon 
gas [l?] ]. The additional term in the chemical potential has the same form as //i introduced 
in the Gibbs relation. This result could be predicted a priori since the physical sense of 
this factor is a variation of an energy of the one liquid particle being due to the phonon 
contribution under changing a number of the particles. 
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Taking into account the relations (|T%|). (J33J), P5|l. and (jSSjl . two flows, // and 9* say, may 
be expressed in terms of n % and S l . Comparing them with (J2J) one can obtain the expressions 
for the coefficients A, B, and C: 

c 2 /c 2 - 1 

A = -fxT -^—i 7> 

fj,n\\ - s\\T (c 2 /c 2 - 1) 

B = n 2 { m -s\\T -I)}' 1 , (55) 



= Tc 2 {^|| + g ||T(l-c^/c 2 )} 

" S ||C 2 {^|| ( C 2/ C 2_ I)}" 



Since superfluid dynamics implies that A < and 5, C > 



have to be positive. These coefficients were estimated in the work 



re denominator in (jB3|) 
in which the phonon 
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gas in superfluids is considered from the statistical standpoint. This result can be derived 



from (|55j) under the condition /my ^> s\\T(c 2 /c 2 ) 



V. SOUND IN SUPERFLUIDS AT LOW TEMPERATURE 
A. Minkowski metrics 

Let us consider sound propagation in superfluids being in a flat space-time at the tem- 
perature close to zero. In this temperature range a free path of the quasi-particle increases 
rapidly and equilibrium in the phonon gas has no time to be established. Under these con- 
ditions, the hydrodynamic equations for the phonon gas are inapplicable. Properties of this 
gas are described by the kinetic equation, in which the collision integral is negligibly small. 
Thus the phonon distribution function satisfies the Liouville equation (|19|) with the right- 
hand side being equal to zero. The Liouville equation must be supplemented by equations of 
motion of the medium where the phonons are situated. These equations are the conservation 
law of the particle number flow (|36|) and the irrotational condition for the superfluid velocity 
©• 

Solution will be sought in the form 

f = fo + fie lhx \ (56a) 
n|| =n + nie iklX \ (56b) 
P' \;/ • r V ' 7 ' ' '". (56c) 

where fi, ni, w J are the small addition to the constant equilibrium values fo, n Q , Vq . 
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The equilibrium distribution function must turn the collision integral into zero. The 
solution of the equation J(x,p) = for the phonon gas (that obeys the Bose-Einstein 
statistics) is well-known (see e.g. Ref. 0) 

/o = eP^-i' (57) 

where f3 k = f3V k } (3 = c(k B T)-\ 

After linearization, the following system of equations is derived from (fTTIJ) . (J3H|) and (jHJ) 

^ Wi + + -te kiV ° { U " J m ^^r ni i = ' (58a) 

(c 2 /c 2 -i)c r 

fclim + n %^ + — y e%7r J /i dP = 0, (58b) 

i <9 /w<9c p \ c 3 dc p f 2 

= + ^ u^j ini y /x ' ( } 

where k\\ = kiVg, Kj = A^ki. In the equation (|58b|) we have taken into account the expression 
(J53J) and the fact that Qi = in equilibrium. 

The equations ()58c|) allow to conclude that Vj is parallel to and only one equation 
remains by entering a new variable v so that Vj = VKj. 

Let us express fi from the equation (|58a|) and substitute it in the integrals in (I58bj) and 
(I58cj) . It is convenient to perform the integration in the superfluids local rest frame, marked 
by the condition V 1 = 5q. In this frame 



9 f^ 2 j ( , c , \ c , _ fc„ <9(c 2 /c-, I I ' 



3) W . ^ -1 



/i = — < [ k\\ k cos $ — kv cos?? — > < — k\\ — k cos?? > , (59) 

oe p [\ c p J c p 2 dn\\ J [c p J 

where fc = (KjK^) 1 ' 2 , p = (iXiix 1 ) 1 ! 2 , and $ is the angle between the vectors Kj and %j. Having 
integrated over the angles we will deduce the system of equations 



<r nil ani| \c p on||/ c re \CpOn\\ ) \ | c/c vcr 




where 



2ttc In Cfc|1 + / p4 dp - 71-2 ^ gT ^ ° 2 In cfc|1 + ° pk (61) 

cfc|| — Cpk J (27rh) 3 de 15 (hc p ) 3 c 2 crey — c p k' 
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Only logarithmic terms are retained after integrating over the angles. They are large in 
comparison with another terms b as we will see below. 

Having equated the determinant of the system ()6U|) to zero we will obtain the dispersion 
relation 



T-2-2 = 1 + — -P L + 1 -^ • ( 62 ) 



The relativistic addition to the sound speed 



Sc p I ( n\\ dc p c p 



c p 2pn\\ { c p on- n2 
differs from the nonrelativistic one 



2 1 2 

+ 1-J (63) 



K = (—X (-^ + iV ^r^^ (64) 



30h 3 p \ Cp J \c p dp J lu — /cc p ' 

by the ratio of the light speed square to the sound speed one in the braces and the production 
pn\\ instead of nonrelativistic mass density p. 

The relativistic and nonrelativistic chemical potentials are related by u = mc 2 + p NR (see 
e.g. Ref. 28(). Thus the classical formula ()64j) takes into account only rest energy of the 
superfluids particles whereas ()63|) includes also thermodynamic contribution in energy. 

To reduce (J63|) to (}64|) it is necessary that the conditions c p C c and p NR mc 2 are 
performed. 

To compare and in ultrarelativistic limit is incorrect in general, since the latter 
expression is pure nonrelativistic. However it may be useful to show how the sound speed is 
changed when nonrelativistic matter transforms to ultrarelativistic one. For the latter the 
sound speed c p = cj a/3 and 

(5c p ) NR 9un||' 

This ratio demonstrates decrease of sound speed addition when relativistic effects become 
essential. 



B. Robert son- Walker metrics 

In this section we consider the expanding universe with Robert son- Walker metrics that 
comprises, at least in part, superfluid matter. Similar models are regarded earlier [29( for 
a Bose-Einstein condensate whose properties intimately related to the superfluidity and 



17 



superconductivity. In particular it is found that the coherence length of the cosmic Bose- 
Einstein condensate is equal to the Jeans wavelength. Thereby the condensate can form 
gravitationally stable structures. If a superfluid state is realized in the cosmic matter, then 
the question naturally arises as to whether the excitations have an influence on this scale. 
Since the Jeans scale varies directly as the speed of sound, evaluation of 5c p may be useful 
to estimate the extent of this influence. 

For the present moment we leave aside the questions of the universe evolution under the 
superfluids influence and here we study behavior of the sound speed near equilibrium. 

Before to study the weakly non-equilibrium state of the phonon gas it is necessary to 
examine a solution of the equations ()19j). (|36|) and Q at equilibrium. 

The phonon gas is supposed to be at rest with respect to the superfluids frame and the 
distribution function depends on only one variable 

fo(x,p) = fo(PkP k ), (66) 

and has the form ()57|) turning the collision integral into zero. For the Robertson- Walker 
metrics with the scale factor a the dispersion relation (J7j) takes the form 

c a 

The kinetic equation leads to a law connecting the variations of the temperature and 
the sound speed in time 

!_^ + ^ = const. (67) 

p a c p a 

The equations of motion (jHo^l and Q of the liquid take the form 

AjVfcfi = fiDVj = 0, (68a) 

DriQ + n 6 = Dn H — -no = 0, (68b) 
a 6 



where the dot denotes the time derivative. It follows from the equations ()68)1 that the 
chemical potential is an arbitrary function of time and the equilibrium particle number 
density of the liquid varies according the law noa 3 =const. 

Let us look for the weakly non-equilibrium solution in the form ()56|). However, unlike the 
flat metrics case, both the equilibrium values f , n o, Vq anc ^ ^ ne sma U additions fi, rii, f- 7 
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depend on time. These quantities vary slower than the exponent power, i.e. 

-7T- < kji, 
ox 1 

and so on. 

Having repeated the procedures of the previous section, we will obtain that 



f 1= Part— <[ku /ccostf — fcvcosfl--^ V w }{ — h\ - - cost? } , (69) 

V I \ ac p J ac p 2 on\\ J I c p a 

and integration over the angles gives the following result 

2n(3 ack\\ + Cpk f p 4 



a 5 ack\\ — Cpk J (2ttH) 
tt 2 (k B T) 4 c 2 ack\\ + c p fc 
15 (hep) 3 c 2 ack\\— Cpk 1 

where the prime denotes the derivative with respect to the argument. To write down the 
dispersion relation for the Robertson- Walker metrics, one should add a 2 to the numerator 
of the right hand side of the expression (1621) . 

In view of (JBTj) and (|68b|) 5c p oc l//xa (logarithmic dependence is neglected, since it related 
with dispersion of phonons). Further consideration is caused by behavior of the chemical 
potential, that can not be determined without additional assumptions. As possible models 
one can consider the simple equations of state that allows to examine the evolution of Sc p . 

(i) Linear barotropic equation of state takes the form P = (c 2 /c 2 )fion\\ for superffuid 
matter at zero temperature where c p = const. Integrating the relation (0J) one obtains 

/c^ o 2 / 2 

fiQ oc n,F and therefore 5c p oc a p' c . 

In the ultrarelativistic limiting case the addition to the sound speed is independent from 
the scale factor a. Thus the sound speed is constant but does not coincide with one at 
zero temperature. For the nonrelativistic limiting case 5c p oc a~ l and universe inflationary 
implies decay of the phonon contribution into the sound speed. 

(ii) Polytropic equation of state P = KnJ^ (7 > 1 is supposed in judgements) in combi- 
nation with the relations (jlj) permits to obtain the expressions: 



Cp 

C 2 1 fc^,,^7-l 



mc 2 + K^nJ, /(7 — 1) 

/i = mc H -An , 

7 — 1 11 
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which is similar to the expressions for nonrelativistic polytropes [3Cj because of their tem- 
perature independence. 

In the nonrelativistic area (mc 2 3> KnJ' 1 ) the sound speed c p oc a -3 ^ -1 -*/ 2 and 5c p oc a -1 . 
One can see that for 7 > 5/3 the sound speed decreases slower then the phonon addition. In 
general, the same behavior also takes place in the relativistic area (mc 2 ~ KnJ' 1 ) especially 
for the large values of a. When mc 2 <C i^njj -1 and the scale factor is not so large, then 
c p = const and 5c p oc a 37 " 4 . The value 7 = 4/3 leads us again to the ultrarelativistic limits 
described in (i). 



VI. CONCLUSION 



In this paper we introduced the principle of constructing macroscopic flows of the phonon 
gas in superfluids using the phonon distribution function. These definitions can be applied 
to arbitrary quasi-particles without any changes. Nevertheless we restrict ourselves by con- 
sideration of the phonon gas. In this case there is no dependence on the quasi-particle energy 
in ifi so this tensor can be used as an effective second metrics. 

This allowed us to apply the well developed technique to derive the macroscopic balance 
equations for the phonon gas, which is the part of the hydrodynamic conservation laws. 
Based on this statement the phonon contribution in superfluid dynamics was obtained. 

It is worth noting that using the expression (j!6b|) instead of (|llb|) . one can draw an 



analogy with the works |3jJ, |32j, |33( (see also references therein) where a distinction between 
momentum and pseudomomentum is discussed. This problem is considered for dielectric 
media and electromagnetic waves in j^Hs^l and for acoustic waves in 33 1. In these papers 
the mentioned problem is studied for an individual optical or acoustic vibrations. In contrast 
to such approach, the tensors (jllbj) and (j!6b|) do not contain parameters of individual 
phonons but describe the quasi-particle ensemble. Therefore, it is reasonable to consider 
concepts of momentum and pseudomomentum in order to compare their applicability to 
description of macroscopic objects. 

Obtained in this paper kinetic equation was solved directly to determine the addition to 
the speed of sound at temperature close to zero. 

The result for the flat space-time is the relativistic generalization of the classical Khalat- 

n 

nikov's one |l5j, whereas results from Sec. IV Bl should be reviewed from the other standpoint. 
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They are certainly preliminary and estimating ones, however they demonstrate that includ- 
ing phonons into cosmological models can be used to fit necessary parameters properly. 

The present research was carried out for the low temperature quasi-particle gas, when 
purely phonon processes dominate, whereas contributions of quasi-particles corresponding to 
another parts of the spectrum are neglected. To develop the similar theory for another sorts 
of quasi-particles, such as rotons, special techniques are needed because of the nonlinear 
dependence between the energy and momentum of a quasi-particle. 
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